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ABSTRACT 

We use highly spectroscopically complete observations of the radio sources from the VLA 1.4 GHz 
survey of the Hubble Deep Field- North region to study the faint radio galaxy population and its evolution. 
The fraction of radio sources that can be optically spectroscopically identified is fairly independent of 
radio flux, with about 60 — 80% identified at all fluxes. We spectrally classify the sources into four spectral 
types: absorbers, star formers, Seyfert galaxies, and broad-line AGNs, and we analyze their properties 
by type. We supplement the spectroscopic redshifts with photometric redshifts measured from the rest- 
frame ultraviolet to mid-infrared spectral energy distributions. Using deep X-ray observations of the 
field, we do not confirm the existence of an X-ray-radio correlation for star-forming galaxies. We also 
do not observe any correlations between 1.4 GHz flux and R magnitude or redshift. We find that the 
radio powers of the host galaxies rise dramatically with increasing redshift, while the optical properties 
of the host galaxies show at most small changes. Assuming that the locally determined far-infrared 
(FIR)-radio correlation holds at high redshifts, we estimate total FIR luminosities for the radio sources. 
We note that the FIR luminosity estimates for any radio-loud AGNs, which we conservatively do not try 
to remove from the sample, will be overestimates. Considering only the radio sources with quasar-like 
bolometric luminosities, we find a maximum ratio of candidate highly-obscured AGNs to X-ray-luminous 
(Lo.5_2 kcv or L2-8 kov > 10'*^ ergs s~-^) sources of about 1.9. Finally, we use source-stacking analyses to 
measure the X-ray surface brightnesses of various X-ray and radio populations. We find the contributions 
to the 4 — 8 keV light from our candidate highly-obscured AGNs to be very small, and hence these sources 
are unable to account for the light that has been suggested may be missing at these energies. 

Subject headings: cosmology: observations — galaxies: active — galaxies: distances and redshifts — 
galaxies: evolution — galaxies: formation 



1. INTRODUCTION 

In the local universe, Seyfert 2 galaxies outnumber 
Seyfert 1 galaxies by 4 to 1 (Maiolino & Rieke 1995), 
and more than half of Seyfert 2 nuclei are Compton- 
thick with obscuring column densities Nh > 10^^ cm^^ 
(Maiolino et al. 1998; Risalti et al. 1999; Guainazzi et al. 
2005; Comastri 2004). In the high-redshift universe, X- 
ray background population synthesis models (refined af- 
ter the discovery of a lower redshift distribution than ex- 
pected for the hard X-ray sources; e.g., Barger et al. 2002, 
2003; Szokoly et al. 2004) similarly predict a population 
of highly-obscured active galactic nuclei (AGNs) that are 
missed in the deep Chandra and XMM-Newton hard X-ray 
surveys (e.g., Franceschini et al. 2002; Gandhi & Fabian 
2003; Ueda et al. 2003; Comastri 2004; GiUi 2004; Fabian 
& Worsley 2004; Treister & Urry 2005; Ballantyne et al. 
2006). Such theoretical predictions appear to be receiving 
some observational confirmation in the work of Worsley et 
al. (2005). Using a source-stacking analysis on the Chan- 
dra Deep Field-North (CDF-N) and Chandra Deep Field- 
South (CDF-S) fields, Worsley et al. (2005) found that the 
resolved fraction in the 6 — 8 keV band was only about 60%. 
(Worsley et al. 2004 did a similar analysis on the XMM- 
Newton Lockman Hole field and found a 50% resolved frac- 



tion above 8 keV.) According to these authors, the missing 
X-ray background (XRB) component has a spectral shape 
consistent with a population of highly-obscured AGNs at 
redshifts z 0.5 — 1.5 having obscuring column densities 
oiNn ~ 10^3 - 1024 cm-2. 

Of course, to derive these results, Worsley et al. (2005) 
had to choose which estimate of the 2 — 8 keV extragalac- 
tic XRB spectrum to use, and there is a great deal of un- 
certainty in these measurements (see Figure 15 of Hickox 
& Markevitch 2006 for a summary; Worsley et al. 2005 
used the XMM-Newton measurement made by Deluca & 
Molendi 2004, but Revnivtsev et al. 2005 found that the 
intensity of the XRB measured by focusing telescopes is 
higher than that measured by coUimated experiments by 
about 10— 15%, possibly due to the extreme complexity of 
measuring effective solid angles of focusing telescopes). In- 
deed, Barger et al. (2002) found that the issue of whether 
there is a need for a substantial population of as-yet un- 
detected highly-obscured AGNs depends critically on how 
the low-energy and high-energy XRB measurements tie 
together. 

Regardless, the combination of the theoretical modelling 
predictions for a population of highly-obscured AGNs 
missing from the deep X-ray surveys and the possible ob- 
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scrvational evidence for a substantial unresolved fraction 
of the XRB has inspired observers to try to quantify the 
fraction of obscured AGNs that may be being missed. 
Deep mid-infrared (MIR) and radio images are an obvi- 
ous avenue for searching for highly-obscured AGNs, since 
extinction in the MIR and radio is small. Various ap- 
proaches using MIR or X-ray data (e.g., Alonso-Herrero 
et al. 2006; Polletta et al. 2006), the combination of MIR 
and radio data (e.g., Martmez-Sansigre et al. 2005), or the 
combination of MIR, radio, and X-ray data (e.g., Donley et 
al. 2005) have been adopted, and a candidate population 
of highly-obscured AGNs has been identified. However, 
the various selection effects in these approaches mean that 
a reliable upper limit on the number of highly-obscured 
AGNs that arc undetected even in the 2 — 8 keV band has 
not yet been obtained. 

In this paper, we take an alternative approach to look for 
highly-obscured AGNs at high redshifts using a pure mi- 
croJansky radio survey selection. Ultradecp radio surveys 
are very useful for tracing the evolution of dust-obscured 
galaxies and AGNs. In contrast to optical surveys, which 
may omit dusty sources, and submillimeter surveys, which 
currently have low spatial resolution and arc limited to 
small areas for any significant depth, the 1.4 GHz surveys 
do not suffer from extinction, provide subarcsecond posi- 
tional accuracy, and cover large areas. 

Although the Jansky and milliJansky radio source pop- 
ulations arc dominated by powerful AGNs, at the mi- 
croJansky level, the radio source population is increas- 
ingly dominated by star-forming galaxies, which produce 
non- thermal radio continuum at 1.4 GHz through syn- 
chrotron emission from supernova remnants. In local star- 
forming galaxies and radio-quiet AGNs, it is well-known 
that the radio power is tightly correlated with the far- 
infrared (FIR) luminosity (e.g., Helou et al. 1985; Condon 
et al. 1991; Condon 1992), probably as a result of both be- 
ing linearly related to the massive star formation rate. If 
we assume that this correlation continues to hold at high 
redshifts (note that there is some observational support 
that the correlation does extend to 2: > 1 from the work 
of Appleton et al. 2004, though the number of galaxies at 
the higher redshifts is relatively low), then we can use this 
relation to estimate the total FIR luminosity of a galaxy 
from its radio power alone. Indeed, even with the advent 
of the Spitzer Space Telescope, this is still the most robust 
way to estimate a galaxy's total FIR luminosity, due to the 
limited MIPS sensitivities at 70 and 160 /im and the uncer- 
tainties in the conversion from 24 /zm luminosity to total 
FIR luminosity. Although estimating total FIR luminosi- 
ties from 24 jim data has been advocated (e.g., Appleton 
et al. 2004; Marcillac et al. 2006), the conversion depends 
on the assumed shape of the spectral energy distribution 
(SED), which is expected to change as a function of dust 
temperature and metallicity. 

Using their Spitzer Infrared Nearby Galaxies Survey 
(SINGS) data, Dale et al. (2005) fomid a total range in 
rest-frame 8 /xm (which corresponds to observed-frame 
24 yum &t z = 2) to total FIR luminosity ratios of more 
than a factor of 20. When they excluded the metal-poor 
dwarfs, the total range dropped to a factor of 10, which 
they argued should be taken as a minimum systematic 
uncertainty for total FIR luminosities inferred from rest- 
frame 8 /um fluxes alone. They cautioned that if the metal 



abundance of the population cannot be inferred a priori, 
then the uncertainties may be considerably higher. Like- 
wise, they found that the rest-frame 24 /xm to total FIR 
ratio spans a factor of 5 in the SINGS sample. 

In addition to our assumption that the FIR-radio cor- 
relation continues to hold at high redshifts, the other pos- 
sible concern with our procedure of estimating the total 
FIR luminosities of the radio sources using the FIR-radio 
correlation is the presence of radio-loud AGNs. For these 
sources, the radio continuum will be in excess of the level 
expected from star formation, and we will overestimate 
their total FIR luminosities. Fortunately, however, since 
we are only interested in obtaining upper limits on the 
number of highly-obscured AGNs, this is not a major con- 
cern for our present analysis, and we have conservatively 
not made any effort to remove such sources from our sam- 
ple. 

Finally, we note that even in the case of the radio-quiet 
AGNs, we are not able to separate out the AGN contri- 
bution to the FIR luminosity from the star-forming con- 
tribution. However, since local radio-quiet AGNs are em- 
pirically observed to follow the FIR-radio correlation, we 
have an upper limit on the AGN FIR emission, which is 
sufficient for our analysis. 

In this paper, we use the ultradecp 1.4 GHz observations 
of the Hubble Deep Field-North (HDF-N; Richards 2000), 
for which we have obtained very complete spectroscopic 
identifications, to study the properties of the radio sample. 
The HDF-N has a wealth of existing multiwavelength data, 
including the HST ACS Treasury Great Observatories Ori- 
gins Deep Survey-North data (GOODS-N; Giavalisco et al. 
2004) and the Spitzer Space Telescope Legacy GOODS-N 
data (P.I. M. Dickinson). In addition, we have recently 
obtained deep wide-area J— and if— band observations 
using the University of Hawaii 2.2 m telescope (Trouille et 
al. 2006, in preparation), which enable us to assign robust 
photometric redshifts to the spectroscopically unidentified 
sources. 

The HDF-N is also the site of the deepest X-ray im- 
age of the sky, the 2 Ms CDF-N. We use these data to 
examine the issue of whether there exists an X-ray-radio 
correlation for star-forming galaxies. In addition, we use 
these data to remove the X-ray-luminous (Xo.5-2 keV or 
-£'2-8 keV > lO'^^ ergs s~^) sources from our radio-identified 
candidate highly-obscured AGN population for tighter up- 
per limits. Finally, we measure the X-ray surface bright- 
nesses for our candidate highly-obscured AGN population 
and others, which we compare with the shape and intensity 
of the XRB. 

The structure of the paper is as follows. In §2, we de- 
scribe the radio sample and the available multiwavelength 
imaging and spectroscopic data. In §3, we discuss the spec- 
troscopic identifications of the radio sample. In §4, we 
briefiy describe our photometric redshift method, which 
uses the optical to MIR SEDs, and we show our redshift 
distribution for the radio sources. In §5, we discuss the 
24 /im properties of the radio sample and the difficulties 
with using these data to estimate FIR luminosities. In §6, 
we describe our spectral classification scheme and exam- 
ine the properties of the radio sample by spectral class. 
In §7, we discuss the X-ray properties of the radio sam- 
ple and investigate whether there is an X-ray-radio corre- 
lation for star-forming galaxies at these higher redshifts, 
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like there is locally. In §8, we describe the optical prop- 
erties of the radio sample. In §9, we determine the FIR 
luminosities of the radio sources using our redshifts and 
the FIR-radio correlation. We then put quantitative up- 
per limits on the numbers of highly-obscured AGNs that 
could have gone undetected in this deepest of hard X-ray 
surveys, and we discuss the implications for the history of 
supermassive black hole growth. In §10, we measure the 
X-ray surface brightnesses of our radio-identified candi- 
date highly-obscured AGN population, as well as of other 
radio and X-ray populations, and we compare them with 
the shape and intensity of the XRB. We summarize our 
results in §11. 

We assume CIm = 0.3, JIa = 0.7, and Hq = 
70 km s-i Mpc-^ All magnitudes are in the AB mag- 
nitude system. 
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2. THE DATA 
2.1. Radio and X-ray Data 

Richards (2000) presented a catalog of 1.4 GHz sources 
detected in the Very Large Array (VLA)^ map of the HDF- 
N, which covers a 40' diameter region with an effective 
resolution of 1.8". The absolute radio positions are known 
to 0.1" - 0.2" rms. 

Alexander et al. (2003) presented the 2 Ms X-ray 
image of the CDF-N, which they aligned with the 
Richards (2000) radio image. Near the aim point, 
the X-ray data reach limiting fluxes of /2-8 kcV ~ 
1.4 X 10"-^^ ergs cm~^ s~-^ and /0.5-2 koV ~ 1-5 x 
10~^^ ergs cm~^ s~^. We determined the maximal radius 
radio circle that fits within the Chandra area to be a 10' ra- 
dius circle centered on the position R.A. = 12''36™54.29*, 
decl. = 62°14'16.0" (J2000.0). This also roughly matches 
the highest sensitivity region of the VLA observations, pro- 
ducing a relatively uniform radio map. Our radio sample 
consists of the 207 radio sources contained within this ra- 
dius. 

For sources within our 10' radius, matching X-ray coun- 
terparts from the Alexander et al. (2003) catalogs to the 
radio sources is not critically dependent on the choice of 
match radius. This can be seen from Figure 7 of Alexander 
et al. (2003) , which shows the positional offset between the 
X-ray and radio sources versus off-axis angle. We chose to 
use a 1.5" search radius as a reasonable compromise be- 
tween not pushing too hard on the accuracy of the data 
and not introducing too much random error, but only two 
additional sources would have been matched with a 2" 
search radius, and only three of the current sources would 
not have been matched with a 1" search radius. 

Richards (2000) gave a 5(t completeness limit of 40 /iJy 
for compact sources in the central region of the map. The 
completeness limit is somewhat higher for extended ra- 
dio sources. We independently determined the complete- 
ness of the Richards (2000) 1.4 GHz catalog by measur- 
ing the radio fluxes of a large, optically-sekxrted sample in 
3" diameter apertures and adjusting the normalization to 
match the measured fluxes in the Richards (2000) cata- 
log for the overlapping sources. Figure la shows that the 
Richards (2000) catalog is highly complete to 70 /xJy and 
then drops to about 50% completeness at 60 /iJy. 

^The VLA is a facility of the National Radio Astronomy Observatory 
operated under cooperative agreement by Associated Universities, Inc. 
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Fig. 1. — Completeness of the HDF-N radio catalog from Richards 
(2000) as a function of 1.4 GIfz flux based on (a) measuring the ra- 
dio fluxes of an optically-selected sample, and (b) comparing the 
Richards (2000) catalog with the Biggs & Ivison (2006) catalog. 
The solid lines show our adopted 60 ynJy completeness limit. 

Biggs & Ivison (2006) have recently given a new cata- 
log of 1.4 GHz sources using all of the A-array data from 
the VLA. This catalog is slightly more sensitive than the 
Richards (2000) catalog and may also be used to estimate 
the completeness of the Richards (2000) catalog. Within 
our 10' radius region, the new catalog contains 280 sources, 
of which 259 have fluxes greater than 40 /iJy, as opposed to 
the 207 sources in the Richards (2000) catalog. Of the 207 
Richards (2000) sources, 188 have direct counterparts in 
the Biggs & Ivison (2006) catalog. For two other sources, 
the positions in the two catalogs are offset due to different 
centerings on the radio sources. The fluxes of the overlap- 
ping sotirces are in reasonable agreement. Most of the 17 
Richards (2000) sources that are missing in the Biggs & 
Ivison (2006) catalog have optical counterparts and red- 
shifts, so they are probably real and have just dropped 
below the flux limit of the Biggs & Ivison (2006) catalog. 
In Figure lb, we show the completeness of the Richards 
(2000) catalog based on the Biggs & Ivison (2006) cata- 
log. The plot shows the ratio of the number of sources 
per flux bin detected in the Richards (2000) catalog to the 
number of sources per flux bin detected in the Biggs & 

(NRAO). The NRAO is a facility of the National Science Foundation 
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Ivison (2006) catalog. Consistent with our analysis of the 
completcincss of the Richards (2000) catalog using the op- 
tical data, this shows that the Richards (2000) sample is 
substantially complete above 70 /xJy but becomes progres- 
sively more incomplete at lower fluxes. In our subsequent 
analysis, we adopt 60 /xJy as our effective completeness 
limit. 

2.2. Optical and Near-Infrared Data 

Giavahsco et al. (2004) presented HST ACS images of 
the HDF-N region (the GOODS-N), and Capak et al. 
(2004) presented ground-based deep optical imaging of a 
very wide-field region encompassing the GOODS-N. The 
ground-based data cover the whole MIPS and IRAC areas 
(§2.3) in the U, B, V, I, z' , and HK' bands. AU of 
the images were registered to the Richards (2000) radio 
catalog for accurate absolute astrometry. 

Trouille et al. (2006, in preparation) carried out deep 
J— and iJ— band imaging of the entire GOODS-N region 
using the Ultra-Low Background Camera (ULBCAM) on 
the University of Hawaii 2.2 m telescope during 2004 and 
2005. ULBCAM consists of four 2kx2k HAWAII-2RG ar- 
rays (Loose et al. 2003) with a total 16' x 16' field of view. 
The images were taken using a 13-point dither pattern 
with ±30" and ±60" dither steps in order to cover the 
chip gaps. The near- infrared (NIR) data were flattened 
using median sky flats from each dither pattern. The im- 
age distortion was corrected using the astrometry in the 
USNO-Bl.O catalog (Monet et al. 2003). The flattened, 
sky-subtracted, and warped images (with typical seeing 
0.7") were combined to form the final mosaic, which has a 
20' X 20' area that fully covers the GOODS-N region. The 
final image was registered to the Richards (2000) radio 
catalog. The integration times at each pixel are 9 hours in 
J and 12.5 hours in H, respectively, and the 5f7 sensitiv- 
ities are 0.84 /iJy and 2.06 /iJy, corresponding to 5(7 AB 
magnitude limits of 24.1 and 23.1, respectively. 

We measured all of the optical and NIR magnitudes at 
the positions of the radio sources using 3" diameter aper- 
tures and corrected them to approximate total magnitudes 
using an average off'set for each waveband (Cowie et al. 
1994). 

2.3. Mid-Infrared Data 

Wang et al. (2006) combined the reduced DRl and 
DR2 IRAC superdeep images from the GOODS-N Spitzer 
Legacy Science Program first, interim, and second data re- 
lease products (DRl, DR1+, DR2; Dickinson et al. 2006, 
in preparation), weighted by exposure time, to form 3.6, 
4.5, 5.8, and 8.0 /im images that fully cover the GOODS- 
N area. We measured the source fluxes at the radio po- 
sitions using fixed 4'.'8 (3.6 and 4.5 /xm) and 6" (5.8 and 
8.0 /xm) diameter apertures. These apertures are approx- 
imately three times the ^ 1'.'7 (3.6 /im) to ^ 2" (8.0 /im) 
FWHM of the point spread function (PSF) and are a good 
compromise between the PSF size and the source sepa- 
ration. We applied aperture corrections from the IRAC 
in-flight PSFs (January 2004), which are consistent with 
the ones published in the IRAC Data Handbook, to the 
measured fluxes. The corrected IRAC fluxes should be 
reasonably close to the total fluxes of the sources, because 
the majority of the sources are point-like compared to the 
~ 2" IRAC PSF. The primary errors in the photometry 



are caused by the high density of sources, especially at 
3.G and 4.5 fim. In these two bands, the typical distance 
between sources is comparable to the PSF, and the maps 
are confusion-limited. Consequently, both the background 
estimate and the aperture photometry are highly subject 
to blending with nearby sources. 

Also following Wang et al. (2006), we directly used the 
DR1+ MIPS 24 /xm source list and version 0.36 MIPS 
24 /im map provided by the Spitzer Legacy Program. This 
source catalog is fliix-limited at 80 /iJy and is a subset of a 
more extensive catalog (Chary et al. 2006, in preparation). 
With the 0.065 deg^ area coverage, the catalog contains 
1199 24 /im sources and is highly complete at 80 /xJy. The 
source positions are based on sources detected in the deep 
IRAC images, and the fluxes are derived using PSF fitting. 
A — 0'.'38 offset in declination was applied to the source 
positions to match the radio-frame astrometry (Richards 
2000). We cross-identified the MIPS sources with the ra- 
dio sources by searching a 1.5" radius region around the 
radio positions. 

2.4. Spectroscopic Data 

Cowie et al. (2004b) and Wirth et al. (2004) report the 
results of an extensive spectroscopic survey of galaxies in 
the ACS GOODS-N region. We supplemented these data 
with spectral observations made with the Deep Extra- 
galactic Imaging Multi-Object Spectrograph (DEIMOS; 
Faber et al. 2003) on the Keck II 10 m telescope subse- 
quent to the publication of these papers and with targeted 
observations of radio and X-ray sources (Barger et al. 2002, 
2003, 2005; Cowie et al. 2004a) made with either DEIMOS 
or the Low-Resolution Imaging Spectrograph (LRIS; Oke 
et al. 1995) on the Keck I telescope. We also use redshifts 
obtained by Chapman et al. (2004, 2005) and Swinbank et 
al. (2004), adopting the Swinbank et al. (2004) NIR red- 
shifts over the Chapman et al. (2004, 2005) optical red- 
shifts, where available, because redshift measurements are 
generally more reliable when they are made from emission- 
line features. Where there are redshifts from both sources, 
the redshifts are within < 0.008 of each other. 

3. RADIO SAMPLE IDENTIFICATIONS 

Cross-identification of the radio sources with their coun- 
terparts is non-trivial. Some radio sources have complex 
structures, and it is only by inspection of overlays of the 
radio contours on optical through MIR images that the 
likely counterparts can be identified. In other cases, the 
counterparts are extremely red, only emerging in the NIR 
or MIR bands. (Reassuringly, these red cores inevitably 
agree with the radio positions.) The most extreme of these 
sources can only be identified in the Spitzer IRAC bands, 
which do not cover the entire radio field. 

There are also cases where previous spectroscopic iden- 
tifications (found through a search of the literature via the 
NASA/IP AC Extragalactic Database, NED) appear to be 
of neighboring, usually much bluer galax;ies, rather than 
of the genuine counterparts. It is possible that some of 
these identifications may, in fact, be correct if the bluer 
structures are parts of a more extended galaxy centered 
on the red core. However, it is very hard to make a con- 
vincing argument that this is indeed the case, rather than, 
for example, that the blue source is a background galaxy 
lensed by the foreground red source. 
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Here wc try to take a conservative approach to the iden- 
tifications. We first determined which sources have un- 
ambiguous counterparts. These include all of the sources 
which have either a 5a ground-based z'- or H-hand coun- 
terpart or an IRAC 5.4 /im counterpart lying within 1" 
of the radio position. This immediately identifies 195 of 
the sources in the sample. The dispersion of these sources 
around the radio positions is 0.3". 

We then overlaid the radio contours of the remaining 12 
sources on the optical images of the field (Fig. 2). One 
source, with a flux of 5.96 mJy, is clearly bi-lobal and cen- 
tered on a small red galaxy. We take this red galaxy to 
be the counterpart to the radio source. The remaining 
11 sources have no obvious identifications, though three 
lie near bright galaxies (one of these could be an off-axis 
source; see below), and a fourth may be associated with 
a pair of interacting galaxies at ^; = 0.642. Nearly all 
of these are faint sources (40 — 80 /xJy), and some may 
be false. (Note that Richards 2000 found a — 9cr source 
within the HDF-N radio image, and since the probability 
of finding such a source in the image is much less than 
1%, he concluded that the noise properties of the image 
are not entirely Gaussian.) However, four of the remaining 
11 sources have fluxes above 70 /iJy, and four (including 
only one of the sources with a flux above 70 /iJy) are in the 
deeper Biggs & Ivison (2006) catalog (see Table 1). These 
may be real sources without counterparts from the optical 
to the MIR. For example, they could be very obscured or 
perhaps high-redshift galaxies. 

Fig. 2. — Mosaic of the bi-lobal source and the 11 remaining 
sources in our radio sample without obvious counterparts. The 
sources are shown in the order of Table 1, starting from the bottom- 
left corner and moving to the right, then looping back to the begin- 
ning of the next row and again moving to the right. The individual 
images arc 12.5" on a side. There are three colors in the images: 
blue is set to be the B image, green is set to be the V image, and 
red is set to be the average of / -|- 2', all from the ACS GOODS-N 
data. The contours show the positions of the radio light. 

Fig. 3. — Mosaic of the 12 sources whose redshifts we omitted 
due to the uncertainty of whether these literature identifications are 
directly associated with the radio sources. The sources axe shown 
in the order of Table 2. starting from the bottom- left corner and 
moving to the right, then looping back to the beginning of the next 
row and again moving to the right. The individual images arc 12.5" 
on a side. There are three colors in the images; blue is set to be the 
B image, green is set to be the V image, and red is set to be the 
average of I + z' , all from the ACS GOODS-N data. The contours 
show the positions of the radio light. 

We summarize the radio and optical coordinates of the 
bi-lobal source and the radio coordinates of the 11 sources 
without obvious conterparts in Table 1. Only one of these 
sources has an X-ray or 24 ^m detection. This source lies 
in the outer regions of a bright z = 0.459 galaxy and may 
be either an off-axis source in that galaxy or a faint back- 
ground source that is hidden by the foreground galaxy and 
hence cannot be identified. 

Next, we determined which of the sources with un- 
ambiguous counterparts had spectroscopic redshifts. We 
omitted redshifts for 12 sources with identifications in the 
literature (from NED), where inspection of the images sug- 
gested that the identifications may not be directly associ- 
ated with the radio sources. As we have noted above, some 
of these identifications may be correct, but it appears safer 
to rely on the photometric redshifts for these sources. We 
summarize the omitted sources in Table 2, and we show 



the radio contours overlaid on the optical images for these 
sources in Figure 3. 
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Fig. 4. — Spectroscopic completeness of our radio sample. The 
sources are grouped into the following flux bins: 40-50, 50—60, GO- 
TO, 70-80, 80-90, 90-100, 100-125, 125-150, 150-200, 200-300, 300- 
500. All of these bins contain 10 or more sources. 

Our final compilation contains 143 spectroscopic red- 
shifts, of which 101 are already in the literature, and 42 
are new to this work. In our subsequent analysis, we dis- 
tinguish the 12 redshifts that we use from Chapman et al. 
(2004, 2005) and Swinbank et al. (2004) from our other 
redshifts, referring to them as the CS sources. In Fig- 
ure 4, we show the fraction of radio sources that are spec- 
troscopically identified. We see that the identified fraction 
is fairly independent of radio flux, with about 60-80% of 
the sources identified at all fluxes. 

25 J 1 
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Flc;. 5. — R magnitude distribution for all of the sources in our 
radio sample (open). Shading denotes the spectroseopically identi- 
fied sample, with darker shading for the CS sources. All seventeen 
R > 27 sources arc placed in the final magnitude bin (hatched). 
Sources with i? < 20 suffer from saturation problems and are likely 
to be brighter than measured. 

In Figure 5, we show the R magnitude distribution for 
all of the sources in our radio sample (open). We de- 
note the spectroscopic sample with light shading and the 
CS sources with dark shading. Beyond R — 24, it be- 
comes more difficult to spectroseopically identify the radio 
sources, while at brighter magnitudes, essentially all of the 
sources can be identified. 
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4. PHOTOMETRIC REDSHIFTS 

Wang ct al. (2006), using the method of Perez-Gonzalez 
et al. (2005), built "training-set" SEDs from 1200 galax- 
ies in the MIPS GOODS-N sample with known rcdshifts 
and spectral types. They built seven templates, ranging 
from an elliptical galaxy spectrum to a very blue star- 
forming galaxy spectrum, over the frequency range from 
6 X lO^'' Hz to 4 X 10^''' Hz. They then made least-squares 
fits of their individual source SEDs to these templates to 
determine the photometric redshifts and spectral types for 
the galaxies in their samples. They found that the method 
worked extremely well over a wide range of redshifts and 
only failed for a small number of sources. Their use of tlic^ 
deep J- and H-h&nd magnitudes and the Spitzer data, in 
addition to the optical data, substantially improved the 
high-redshift end of the photometric redshift determina- 
tions. 

To derive photometric rcdshifts for our radio sources, 
we similarly constructed SEDs using the optical data, the 
Spitzer IRAC 3.6 /zm, 4.5 /xm, 5.8 /zm, and 8.0 /xm data, 
the Spitzer MIPS 24 fim data, and the ULBCAM NIR J 
and iJ-band data. We then made least-squares fits to the 
Wang et al. (2006) templates to determine the photomet- 
ric redshifts and spectral types for the galaxies in our radio 
sample that were detected in at least five bands. This gives 
us an additional 47 redshifts. Thus, in total, only 17 of our 
radio sources (including the 11 without any obvious coun- 
terparts) have neither a spectroscopic nor a photometric 
redshift. We show a comparison of the photometrically 
identified sources (open) with the spectroscopically iden- 
tified sources (light shading) in Figure 6. The CS sources 
are denoted by dark shading. The gap in the spectroscopic 
redshift distribution between ^ ~ 1.4 and z ~ 1.9, known 
as the "redshift desert" , refiects the difficulty of identifying 
galaxies with redshifts in this range, where [OH] 3727 A 
has moved out of the optical window and Lya 1216 A has 
not yet entered in. 
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Fig. 6. — Number of spectroscopically (shaded) and photomet- 
rically (open) identified sources vs. redshift for the radio sample. 
The CS sources are denoted by dark shading. 



5. 24 MICRON PROPERTIES OF THE RADIO SAMPLE 

Studies based on the Infrared Space Observatory (ISO) 
15 fj,m data found a loose correlation between the MIR 
emission and the radio continuum (Cohen et al. 2000; El- 



baz et al. 2002; Garrett 2002; Gruppioni et al. 2003). Ap- 
pleton et al. (2004) confirmed the existence of this corre- 
lation with Spitzer data, and Marcillac et al. (2006) pro- 
posed using it to derive FIR luminosities. However, while 
obtaining the FIR luminosity from a radio measurement 
only depends on the FIR-radio correlation and the radio 
spectral index (and the assumption that the source is not 
a radio-loud AGN), obtaining the FIR luminosity from a 
MIR measurement strongly depends on the library of tem- 
plate SEDs that is used to i^-correct the data. 
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Fig. 7. — Ratio of 24 /im flux to 1.4 GHz flux vs. redshift for the 
spectroscopically identified radio sources in the MIPS area. Solid 
(open) symbols denote sources with (without) 24 /xm detections. 
Circles denote the CS sources. Large open squares denote sources 
with either Lo.B-2 kev oi' ^2-8 keV ^ 10*^ ergs s~^. Large open di- 
amonds denote spectroscopically classified absorbers. Spectroscop- 
ically unidentified sources are nominally plotted at z = 0. The 
overplotted curves are the ratios from redshifted template SEDs of 
M82 (solid), Arp 220 (dotted), and a spiral galaxy (dashed) (Silva 
et al. 1998). 

We illustrate the difficulty with using MIR data to esti- 
mate FIR luminosities in Figure 7 (similar to Figure 1 of 
Donley et al. 2005), where we show the observed 24 y^m to 
1.4 GHz flux ratio versus redshift. Solid (open) symbols 
denote sources with (without) 24 /zm detections, and large 
open squares denote sources with soft or hard X-ray lumi- 
nosities in excess of 10^^ ergs s""'^. For the sources without 
24 /xm detections, we have adopted the 80 /xJy flux limit of 
the 24 /xm source catalog (sec §2.3). We have overplotted 
24 /xm to 1.4 GHz flux ratios versus redshift for redshifted 
template SEDs (Silva et al. 1998). It is clear that there is 
substantial spread in the ratio, with the 24 /xm measure- 
ments being very sensitive both to redshift and to galaxy 
type. For the most part, the data points have values that 
are within the range expected from the redshifted tem- 
plate SEDs; however, there are also a number of sources 
not detected at 24 /xm, including many of the absorbers 
{large open diamonds; sec §6). The absorbers are most 
likely radio-loud AGNs, which would not be expected to 
have SEDs like those shown. 

We note that Donley et al. (2005), in their search for 
highly-obscured AGNs, defined a selection threshold of 
q = log(S'24 pm/5'i.4 GHz) < to classify a galaxy as prob- 
ably having an AGN. Based on the SEDs used in our Fig- 
ure 7, this criterion might result in some contamination 
by ultraluminous infrared galaxies, and, to a lesser extent, 
by spirals, as both of those SEDs dip below g = 0. We do 
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sec sources below this value that do not have X-ray lumi- 
nosities typical of an AGN and arc not spcctroscopically 
classified as absorbers, and some of these sources may be 
highly-obscured AGNs, but clearly we cannot rely on the 
24 /xm data to determine their FIR luminosities. 

6. PROPERTIES OF THE RADIO SAMPLE BY SPECTRAL 

CLASS 

For our spectroscopically identified sources, we classi- 
fied the optical spectra into four spectral classes, roughly 
following the procedure used by Sadler et al. (2002) to an- 
alyze low-redshift 1.4 GHz samples. We classified sources 
without any strong emission lines [EW([OII])< 3 A or 
EW(Ha+NII)< 10 A)] as absorbers; sources with strong 
Balmer lines and no broad or high-ionization lines as star 
formers; sources with [NeV] or CIV lines or strong [OIII] 
[EW([OIII] 5007 A) > 3 EW(H/3)] as Seyfert galaxies, 
and, finally, sources with optical lines having FWHM line 
widths > 2000 km s~^ as broad-line AGNs. We have not 
classified one high-redshift source at z = 2.2032. This 
source must be a strong star former because it has UV ab- 
sorption lines; however, its [Oil] emission is redshifted into 
the NIR and hence cannot be measured with our spectra. 
We have also not attempted to spectroscopically classify 
the CS sources, which are generally of this type. 

6.1. Radio-to- Optical Ratios by Spectral Class 

Some radio quasar studies have found bimodal distri- 
butions of rest-frame radio-to-optical ratios (e.g., Kellcr- 
mann et al. 1989; Stocke et al. 1992; Ivezic et al. 2002; 
however, see White et al. 2002; Cirasuolo et al. 2003). 
These have been interpreted as evidence that quasars come 
in two distinct populations, 'radio-loud' and 'radio-quiet'. 
The boundary value between radio-loud and radio-quiet 
has been inferred to be /s ghz/J °a — (Stocke et al. 

2500 j4 

1992). 

For our study, we define a rest-frame optical-to-radio 
ratio of i? = /1.4 ghz/ f^^^Q ^- This is an approximately 
equivalent ratio, since the .ftT-correction to go from 5 GHz 
to 1.4 GHz is very similar to the X-corrcction to go from 
2500 A to 6500 A, assuming spectral indices of a = 0.8 
for fi, oc in both cases (Yun et al. 2001; Zheng et al. 
1997). 

In Figure 8, we show a histogram of logarithmic R for 
our z < 1.6 spcctroscopically classified radio sample. We 
have such uniform and complete wavelength coverage for 
all of the radio sources that we have just interpolated be- 
tween our measurements to obtain rest-frame AB 6500 A 
magnitudes. We have restricted the plot to sources with 
z < 1.6 so that our interpolation remains valid. We have 
also JC-corrected the radio measurements to rest-frame 
1.4 GHz. The star formers clearly dominate the micro- 
.Jansky sample, as was inferred from previous work (e.g., 
Windhorst et al. 1995; Richards et al. 1998; Richards 2000; 
Roche et al. 2002). 

Unfortunately, we do not have a large enough sample 
of Seyfert galaxies plus broad-line AGNs to investigate 
whether there is bimodality in the rest-frame radio-to- 
optical ratios. However, we can see that all of the spec- 
tral classes show a wide range in radio-to-optical ratios. 
Thus, in contrast to the conclusions of Machalski & Con- 
don (1999) from their study of milliJansky radio sources 



in the Las Campanas Redshift Survey, we do not find this 
ratio to be a useful discriminator between star-forming 
galaxies and AGNs. Afonso et al. (2005) came to a similar 
conclusion using the Phoenix Deep Survey, which reaches 
into the sub- 100 /iJy regime. They argued that the radio 
morphology and polarimetry used by Machalski & Condon 
(1999) to identify AGNs results in the misclassification of 
many low-power AGNs, and these sources would blur the 
distinction if properly classified as AGNs. Since the radio- 
to-optical ratio docs not appear to have much diagnostic 
power for our sample, we do not use it subsequently. 
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Fig. 8. — Histogram of the logarithmic rcst-framc radio-to- 
optical ratio, R, for our spcctroscopically classified 2 < 1.6 radio 
sample divided into absorbers, star formers, and Seyfert galaxies 
plus broad-line ACNs (BLAGNs). Each tickmark on the y-axis 
represents one source. The dashed vertical line shows the classical 
boundary (R = 10) between radio-quiet and radio-loud. 

6.2. Radio Power and Redshifts by Spectral Class 

In Figure 9, we show histograms of the (a) logarithmic 
radio powers and (b) redshifts for the spectroscopically 
classified radio sample (shaded). We calculate the rest- 
frame radio powers from the equation 

P1.4 GHz - 4^rfL'5i.4 GHzlO-'"(l + z)"'' crgs s'^ Hz'^ . 

(1) 

Here dz, is the luminosity distance (cm), 51.4 qhz is the 
1.4 GHz flux density (/uJy), and a is the radio spectral 

index, which we take to be 0.8 (Yim et al. 2001). 

The bulk of the star formers tend to lie at the lower radio 
powers. The redshift distribution for the Seyfert galaxies 
plus broad-line AGNs stretches to higher redshifts than 
the star formers and the absorbers, though the absorbers 
in particular become more difficult to identify spectroscop- 
ically at z > 1. We show that this is a selection effect by 
including on Figure 9b the redshift distribution for the 
photometrically identified absorbers classified using their 
SEDs (see §4) '(open). 

7. X-RAY PROPERTIES OF THE RADIO SAMPLE 

An X-ray-radio correlation for star-forming galaxies has 
been determined locally by a number of authors (e.g., 
Bauer et al. 2002; Ranalli et al. 2003; Grimm, Gilfanov, 
& Sunyaev 2003; Gilfanov, Grimm, & Sunyaev 2004), sug- 
gesting that the X-ray emission in star-forming galaxies 
can be used as an indicator of the star formation rate. 
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These same authors, using a smaU sample of sources in 
the CDF-N (at the time, only the 1 Ms Chandra data were 
available; Brandt et al. 2001), also claim that the linear 
relation extends to higher redshifts {z < 1.3). With our 
highly spectroscopically complete, cleanly selected sample 
from just one flux band (1.4 GHz), we are in an excellent 
position to test whether the local X-ray-radio correlation 
does indeed continue to hold to higher redshifts. 
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Fig. 9. — (a) Logarithmic power and (b) rcdshift distributions for 
our spectroscopically classified radio sample (shaded) divided into 
absorbers, star formers, and Seyfert galaxies plus broad-line AGNs 
(BLAGNs). Each tickmark on the y-axis represents one source. In 
(b), we also show the rcdshift distribution for photometrically iden- 
tified absorbers classified using their SEDs. (sec §4) (open). 

In Figure 10a (10b), using different symbols for the dif- 
ferent spectral classes, we plot 2 — 8 keV (0.5 — 2 keV) lu- 
minosity over 1.4 GHz power (we use the ratio to remove 
the redshift effect and improve the dynamic range) versus 
1.4 GHz power for the spectroscopically identified galax- 
ies in our radio sample that were significantly detected by 
Alexander et al. (2003) in the 2 - 8 keV (0.5 - 2 keV) 
band. We calculated the rest-frame 1.4 GHz powers using 
Eq. 1, and we calculated the rest-frame hard (soft) X-ray 
luminosities from 



Lx=4Trdifx{l + zy-^ evgss 



(2) 



Here cIl is the luminosity distance (cm), fx is the 
observed- frame hard (soft) X-ray flux (ergs cm~^ s~^), 
and r is the photon index, which we take to be 1.8 for 



all of the sources. Note that using the individual pho- 
ton indices (rather than the universal power-law index 
of r = 1.8) to calculate the K-corrections would result 
in only a small difference in the rest-frame luminosities 
(Barger et al. 2002). 

We assume that any source more X-ray luminous than 
10^^ ergs s^^ is very likely to be an AGN on energetic 
grounds (Zezas et al. 1998; Moran et al. 1999), so we illus- 
trate the ratio of this fixed hard (soft) X-ray luminosity 
to the 1.4 GHz power versus the 1.4 GHz power with a 
diagonal line in the figures. We note that since none of 
the star-forming galaxies (solid diamonds) show any ob- 
vious AGN signatures in their optical spectra, those with 
Lx > 10^^ ergs are most likely to be obscured AGNs. 

With a dashed horizontal line, we show the local hard 
(soft) X-ray/radio relation found by Ranalli et al. (2003; 
their Eqs. 13 and 9, respectively). We have converted their 
2 — 10 keV relation to 2 — 8 keV, assuming a photon index 
of r = 1.8. Although in Figure 10a we have too few sig- 
nificantly hard X-ray detected, star-forming galaxies with 
L2-S keV < 10^^ ergs s~^ (large open diamonds) to iden- 
tify any kind of correlation, in Figure 10b, we can see the 
correlation that previous authors had identified. However, 
several radio-powerful star formers significantly detected 
in soft X-rays lie well below the apparent correlation and 
may be the luminous edge of the X-ray undetected pop- 
ulation. In fact, it is important to stress that the above 
figures only include X-ray detected sources, not upper lim- 
its. Including the upper limits would be the true test of 
whether a correlation actually exists or whether one is just 
observing the upper end of the distribution and hence be- 
ing fooled into thinking that there is a correlation when, 
in actuality, it is merely a selection effect. 

Unfortunately, from Figures 10c and d, we can see that 
the apparent correlation is indeed merely a selection effect. 
Here we plot only the spectroscopically identified star for- 
mers that are not significantly detected in the 2 — 8 keV 
and 0.5 — 2 keV bands, respectively, using the hard and 
soft X-ray detection limits from Alexander et al. (2003) 
(downward-pointing arrows). The diagonal and horizontal 
lines are as before. The first thing to note is that even 
with the 2 Ms exposure, we do not have the sensitivity to 
measure the bulk of the star-forming galaxy population in 
X-rays. The second thing to note is that — again seen most 
clearly in the soft X-ray figure — the upper limits are not 
consistent with the Ranalli et al. (2003) relation. 

We note that some of the sources in Figure 10 may be 
radio-loud AGNs, but anyone using the supposed corre- 
lation with their X-ray data is going to be faced with 
that uncertainty. Thus, the fact that we do not confirm 
the existence of an X-ray radio correlation at these red- 
shifts, even with our knowledge of the optical spectra of 
the sources, is a major concern for those who want to infer 
star formation rates for their galaxies. Although there are 
attractions to measuring a star formation history based 
on X-ray surveys becatise it avoids problems of obscura- 
tion (e.g., Norman et al. 2004), we are forced to conclude 
that it is not easy to do so. 

8. OPTICAL PROPERTIES OF THE RADIO SAMPLE 

In Figure 11a, we plot 1.4 GHz fiux versus R magnitude 
for the radio sample. Interestingly, there is no observed 
correlation between 1.4 GHz flux and R magnitude, which 
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Fig. 10. — Rest-frame 2 — 8 keV or 0.5 — 2 kcV luminosity over 1.4 GHz power vs. 1.4 GHz power for the spcctroscopically identified 
galaxies in the radio sample (absorbers — solid squaTes; star formers — solid diamonds; Seyfert galaxies — solid upside-down triangles; broad- 
line AGNs — large, solid triangles; CS sources and our unclassified z = 2.2032 source [see §2.4] — solid circles) that are (a) significantly hard 
X-ray detected, (b) significantly soft X-ray detected, (c) star-forming galaxies only without significant hard X-ray detections, (d) star-forming 
galaxies only without significant soft X-ray detections. In (a) and (b), the large, open diamonds denote significantly X-ray detected, star- 
forming galaxies with either (a) L2-8 keV < lO^'^ ergs or (b) io.6-2 keV < 10^^ ergs s~^. In (c) and (d), the arrows show the hard and 
soft X-ray detection limits from Alexander et al. (2003). The dashed horizontal lines show the Ranalli et al. (2003) best-fit relations for their 
local sample assuming a linear slope (their Eqs. 9 and 13 for the soft and hard X-ray bands, respectively; we have K-corrected the latter to 
2 — 8 keV from 2 — 10 keV assuming T = 1.8). The solid diagonal lines assume a fixed X-ray luminosity of 10*^ ergs s~^, the luminosity above 
which the sources are most certainly powered by AGNs. 



10 



helps to explain the uniformity of the spectroscopically 
identified fraction with radio flux (see Fig. 4). We illus- 
trate this in histogram form in Figure lib, where we plot 
number versus R magnitude for all of the radio sources in 
our sample, divided into five radio flux bins. The R mag- 
nitude distribution is consistent with being drawn from 
an invariant population. This result differs from that of 
Georgakakis et al. (1999) and Afonso et al. (2005), who 
found fainter median R magnitudes with decreasing radio 
flux density. The most likely reason for the discrepancy is 
that the optical observations used by those authors were 
much shallower, probing only to ii = 22.5. 
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Fig. 11. — (a) 1.4 GHz flux vs. R magnitude for the radio 
sample (absorbers — solid squares; star formers — solid diamonds; 
Seyfcrt galaxies — solid upside-down triangles; broad-line AGNs — 
large, solid triangles; CS sources and our unclassified z=2.2032 
source [see §2.4] — solid circles; photometric redshifts-open squares; 
unidentified sources — open leftward-pointing triangles). Sources 
with radio fluxes greater than 2500 fiJy are shown at that flux, 
and sources with R > 27 are shown at that magnitude. Sources 
with R < 20 suffer from saturation problems and are likely to be 
brighter than measured. Radio flux ranges in microJanskys are 
given for each bin. Each tickmark on the y-axis represents five 
sources. Darker shading denotes sources with R > 27. 



In Figure 12a, we show radio flux versus redshift for 
the spectroscopically (solid squares) and photometrically 
(open squares) identified sources in the radio sample. We 
plot the unidentified sources at z < 0. Again we see no cor- 
relation between the two variables. In Figure 12b we show 
this in histogram form, plotting number versus redshift for 



all of the radio sources with either spectroscopic or photo- 
metric rcxlshifts, divided into five radio flux bins. The red- 
shift distribution is consistent with being invariant with ra- 
dio flux, although the selection of optically bright sources 
through the spectroscopic identification process may mean 
that we are missing the higher redshift tail of the redshift 
distribution function. Thus, probing to fainter radio fluxes 
does not sample to higher redshifts (Condon 1989). 
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Fig. 12. (a) 1.4 GHz flux vs. redshift for the radio 
sample (absorbers — solid squares; star formers — solid diamonds; 
Seyfert galaxies — solid upside-down triangles; broad-line AGNs — 
large, solid triangles; CS sources and our unclassifled z=2.2032 
source [see §2.4] — solid circles; photometric redshifts-open squares; 
unidentifled sources — open leftward-pointing triangles at z < 0). 
Sources with radio fluxes greater than 2500 ^t.ly are shown at that 
flux, (b) Redshift distribution for the radio sample with spectro- 
scopic or photometric redshifts, divided into five fiux bins. Radio 
flux ranges in microJanskys are given for each bin. Each tickmaxk 
on the y-axis represents five sources. 

This invariance in both the optical apparent magnitudes 

and the redshifts of the host galaxies is a striking result. 
In Figure 13, we plot rest-frame AB 4500 A - 8500 A 
color versus (a) redshift and (b) 1.4 GHz flux for the spec- 
troscopically (solid symbols) and photometrically (open 
squares) identified radio sources. Again, our uniform and 
complete wavelength coverage for all of the radio sources 
means that we can just interpolate between our measure- 
ments to obtain the rest-frame colors. We only use red- 
shifts up to 2 = 1.2 so that our interpolation remains valid. 
The solid and dotted lines show the rest- frame AB 4500 A 
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— 8500 A color of an elliptical and irregular galaxy, respec- 
tively, using Coleman et al. (1980). In both plots, there is 
very little color variation with redshift. The radio sources 
are clearly drawn from the same host galaxy population, 
and the host galaxy properties are not evolving much with 
redshift. 
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Fig. 13.— Rest-frame AB 4500 A - 8500 A color vs. (a) redshift 
and (b) 1.4 GHz flux for the z < 1.2 spcctroscopically (absorbers — 
solid squares; star formers — solid diamonds-. Scyfcrt galaxies — solid 
upside-down triangles; broad-line AGNs — large, solid triangles) and 
photometrically (open squares) identified radio sample. The solid 
(dotted) line shows the rest-frame color of an elliptical (irregular) 
galaxy from Coleman et al. (1980). 



In Figure 14a, we plot the absolute rest-frame R magni- 
tude, Mr, versus redshift for the spectroscopically (solid 
symbols) and photometrically (open squares) identified ra- 
dio sources, interpolating between our measurements to 
obtain rest-frame AB 6500 A, and restricting to 2 < 1.6 
so that our interpolation remains valid. The approximate 
effect of our i? ~ 24 spectroscopic selection limit on Mj^ 
is shown by the dotted curve. The median absolute mag- 
nitude at z < 0.4 is Mji = —21.4, while that between 
z = 0.4 and z = 0.8 is Mr = —21.5, so there may be a 
weak evolution to brighter magnitudes at higher redshift. 
However, this evolution is small, particularly when pas- 
sive evolution of the galaxies is taken into account. At no 
redshift do the host galaxies become much brighter than 
Mr = — 23. Thus, we see a narrow range in Mr with red- 
shift {dashed lines), which suggests that the radio sources 



are chosen from approximately L* optical galaxies at all 
redshifts. At a very crude level, the radio hosts are stan- 
dard candles in their optical properties. 
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Fig. 14.— Absolute AB 6500 A magnitude versus (a) redshift and 
(b) 1.4 GHz power for the z < 1.6 spcctroscopically (absorbers — 
solid squares; star formers — solid diamonds; Seyfert galaxies — solid 
upside-down tria/ngles: broad-line AGNs — large, solid triangles) and 
photomctricalh' {open squares) identified radio sample. Dashed 
lines show the approximate range in absolute magnitude, from —20 
to —23, where the majority of the radio sources lie. Dotted curve in 
(a) illustrates the effect of an il ~ 24 spectroscopic incompleteness 
limit on the absolute magnitudes. 

In Figure 14b, we plot Mr versus 1.4 GHz power. Even 
though the galaxy hosts have similar optical luminosities 
at all redshifts, we can see from the figure that they have a 
wide range of radio powers (see also Cirasuolo et al. 2003). 
Thus, the radio powers of the host galaxies must be ris- 
ing dramatically with increasing redshift, while the optical 
properties of the host galaxies are not changing. 

We can see this dramatic evolution in Figure 15, which 
shows the spcctroscopically (solid symbols) or photomet- 
rically [open squares) identified radio sample in a redshift 
versus 1.4 GHz power plot. The vertical dashed and solid 
lines show the radio powers (see §9 for how we determined 
these) corresponding to a luminous infrared galaxy (LIRG; 
10^"'^ Lq < LpiR < Lq) and an ultraluminous in- 

frared galaxy (ULIRG; LpiR > 10^^ ^o), respectively. 
All of the ULIRGs in the 1.4 GHz sample to a complete- 
ness limit of 60 fxJy (dotted curve) should be detected to 
z ~ 1.3. 
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Interestingly, wc see from Figure 15 that only a narrow 
range of luminosities is observed at each rcdshift and that 
very high-luminosity sources are primarily found at high 
redshifts (2 ^ 1). We therefore conclude that similar opti- 
cal galaxies (see Fig. 14a) are hosting more powerful radio 
galaxies at higher redshifts. (See Cowie et al. 2004a for 
a quantitative description of the evolution of the 1.4 GHz 
luminosity function.) 
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Fig. 15. — Redshift vs. 1.4 GHz luminosity for the spectro- 
scopically (absorbers — solid squares; star formers — solid diamonds; 
Seyfert galaxies — solid upside-down triangles; broad-line AGNs — 
large, solid triangles; CS sources and our unclassified z = 2.2032 
source [see §2.4] — solid circles) and photometrically {open squares) 
identified radio sample. The dashed and solid vertical lines show the 
equivalent radio powers of a LIRG and a ULIRG, respectively. The 
dotted curve shows the radio powers corresponding to the 60 //Jy 
completeness limit of the radio sample. 
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Fig. 16. — 1.4 GHz flux versus 1.4 GHz power for the spectroscop- 
ically identified radio sample, including the CS sources {z < 0.5 — 
solid squares; 0.5 < z < 1 — solid triangles; 1 < 2 < 2 — open 
squares; z > 2 — open circles). The dotted line shows the 60 //Jy 
completeness limit. Sources with radio fluxes greater than 2500 /iJy 
are shown at that flux. 

We expand on this result in Figure 16, where we show 
1.4 GHz flux versus 1.4 GHz power. Here we denote the 

different redshift intervals by different symbols (z < 0.5, 
solid squares; 0.5 < 2; < 1, solid triangles; 1 < z < 2, open 
squares; z >2, open circles). We see a series of increasing 
luminosity slices with increasing redshift. The reason for 
this eff'ect was first suggested by Condon (1989) based on 
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the radio source counts. At low redshifts {z <^ 1), the ra- 
dio surveys are dominated by the most luminous sources 
in each redshift shell because the radio population is ex- 
periencing such a dramatic luminosity evolution. This is 
unlike any other wavelength survey, where sources are gen- 
erally picked up throughout the whole volume distribution. 
However, at higher redshifts, the evolution slows, and a 
wider range of luminosities is observed. 

Figure 16 contradicts the suggestion by Ciliegi et al. 
(2005) that most of the faintest radio sources may be as- 
sociated with relatively low radio luminosity sources at 
relatively modest redshifts rather than with high radio lu- 
minosity AGNs at high redshifts. In fact, the faintest radio 
sources have a wide range of radio powers and redshifts, 
with the higher radio powers corresponding to the higher 
redshifts. 

9. UPPER LIMITS ON HIGHLY-OBSCURED AGNS 

An important goal of the present work is to determine 
an upper limit on the number of AGNs with quasar-like 
bolometric luminosities (~ 4 x 10^^ ergs s~^; see below) 
that could be highly obscured — such that the light emerges 
in the FIR and hence not be detected in current deep X- 
ray surveys. This determination will be an upper limit, 
since for any of the radio sources, it is likely that much of 
the FIR light will be due to star formation rather than to 
AGN activity. Moreover, such a determination will con- 
tain radio-loud sources, where the FIR-radio correlation 
substantially overestimates the FIR. 

We calculate total FIR luminosities for the radio 
sources, assuming that the FIR-radio correlation for star 
formers and radio-quiet AGNs holds at high redshifts, 
and assuming that the sources arc well described by this 
correlation. Following Barger et al. (2001), we use the 
FIR-radio correlation given in Sanders & Mirabel (1996) 
with q = 2.35 that they find holds for sources cover- 
ing several orders of magnitude in FIR luminosity. Then 
fpiR = 8.4 X 10^^ /f 4, where fl^^ is the rest-frame flux at 
1.4 GHz. Assuming a synchrotron spectrum with a spec- 
tral index of 0.8 (Yun et al. 2001), we can calculate the 
total FIR luminosity from the observed 1.4 GHz flux using 



Lfir = 47rdi(8.4 x 10")/i.4(l + z) 



-0.2 



(3) 



where dL is the luminosity distance in cm and /1.4 has 
units of ergs cm^^ s^^ Hz"-'^. 

Of particular interest for our analysis are the sources 
that do not also have high (> 10^^ ergs s~^) X-ray lumi- 
nosities. Quasars are typically defined as having 2 — 8 keV 
X-ray luminosities > 10''^ ergs s""'^, which, with a fac- 
tor of 35 bolometric correction (e.g., Elvis ct al. 1994; 
Kuraszkiewicz et al. 2003), would roughly correspond to 
ULIRG luminosities. Thus, we consider any ULIRGs that 
do not have high X-ray luminosities to be our primary 
candidates to contain highly-obscured AGNs. 

In Figure 17, we plot redshift versus total FIR luminos- 
ity for the sources with ULIRG (vertical line) luminosities 
Lfir > 4 x 10^^ ergs s~^. We denote sources with spec- 
troscopic redshifts by solid symbols and sources with pho- 
tometric redshifts by open squares. We circle the sources 
that have either hard (2-8 keV) or soft (0.5 - 2 keV) X- 
ray luminosities > 10"^^ ergs s~^. The dotted curve shows 
the 60 iiJy completeness limit for the radio observations. 
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The dashed horizontal hne shows an adopted upper red- 
shift cut-ofF of z = 2, which is a compromise to keep out 
ULIRG selection as uniform as possible (since the radio 
completeness limit begins to cut out the lower luminosity 
ULIRGs at z ^ 1.3), while still retaining a statistically 
significant sample of sources. 
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Fig. 17. — Rcdshift vs. total FIR luminosity for the radio sample 
with ULIRG luminosities. Solid squares and solid circles (the latter 
are CS sources) denote spectroscopic rcdshifts, and open squares 
denote photometric rcdshifts. Large, open circles denote sources 
that have either io.5— 2 keV or ^2— 8 keV ^ lO^'^ ergs s^^. The dot- 
ted curve shows the 60 juJy completeness limit of the radio sample. 
The solid vertical line shows the minimum luminosity of a ULIRG. 
The dashed horizontal line shows an upper redshift cut-ofT of 2 = 2 
to provide as uniform a selection of ULIRGs as possible, while still 
retaining a statistically significant sample. 

In Figure 18a. we show histograms of the logarithmic to- 
tal FIR luminosities of the z <2 radio sources with either 
spectroscopic or photometric redshifts and ULIRG lumi- 
nosities divided into two categories: high X-ray luminos- 
ity sources (either Lo.5_2 kcV or L2-8 keV > 10''^ ergs s~^) 
and low X-ray luminosity sources (both Lo.5-2 koV and 
L2-8 keV < 10^^ ergs s~^). There are 38 sources in the low 
X-ray luminosity histogram and 20 sources in the high X- 
ray luminosity histogram. Excluding the bi-lobal source 
(the source with the highest FIR luminosity estimate in 
Fig. 18, which is likely a substantial overestimate for this 
radio-loud AGN) from being considered a candidate to 
contain a highly-obscured AGN, we find a maximum ratio 
of highly-obscured AGN candidates to high- luminosity X- 
ray sources of 1.9 for radio sources in the quasar class. We 
stress again that this is an upper limit, since it includes 
star-forming ULIRGs as well as radio-loud AGNs, where 
the FIR-radio correlation is overestimating the FIR lumi- 
nosity. Even if all 17 sources without photonic;tric; rc;dsliifts 
(11 of which do not have any obvious counterparts) were 
assumed to have ULIRG luminosities, to lie a,t z <2, and 
not to be X-ray luminous, the ratio would only increase to 
2.7. 

To show the classes of sources that are comprising the 
low X-ray luminosity population, in Figure 18b, we divide 
the population into absorbers, star formers, and "other" . 
The latter category includes spectroscopically unclassified 
sources and sources with only photometric redshifts. In 
principle, we should be able to exclude the absorbers, 
which are presumably radio-loud AGNs in the centers of el- 
liptical galaxies, from being considered candidates to con- 



tain highly-obscured AGNs. However, in practice, apart 
from the bi-lobal source, we leave the absorbers in because 
of the difficulties in distinguishing between radio-loud and 
radio-quiet AGNs, as discussed in §6.1. 
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Fig. 18. — (a) Histograms of logarithmic FIR luminosity for the 
LpiR. > 4x 10^'' ergs s^-'^ and z < 2 radio sample with spectroscopic 
or photometric rcdshifts, divided into two X-ray luminosity bins: 
high X-ray (cither L().5_2 kcV ^2—8 koV ^ Iff*^ ergs s^-'-) and low 
X-ray (both Lo.5_2 koV and L2-8 keV < 10'*^ ergs s^-*-). Each tick- 
mark on the y-axis represents five sources, (b) Histograms of the low 
X-ray luminosities from (a) divided into three categories: absorbers, 
star formers, and "other" (spectroscopically unclassified sources and 
sources with only photometric redshifts). The CS sources are de- 
noted by dark shading. Each tickmark on the y-axis represents one 
source. 



In Figure 19, we show in two redshift bins the num- 
ber densities of radio sources with either spectroscopic or 
photometric redshifts and LpjR > 4 x 10^''' ergs s^^, di- 
vided according to X-ray luminosity. At these high FIR 
luminosities, the ratio of highly-obscured AGN candidates 
to high X-ray luminosity sources does not appear to be 
changing much with redshift, with values of 2 and 1.8 for 
the lower and higher rcdshift bins. 

We can get a rough check on whether the upper bound 
on the ratio continues to stay about the same at lower 
FIR luminosities, and hence lower redshifts, by redoing 
our analysis for LIRGs. Of course, there is much more 
contamination at these lower luminosities by star forma- 
tion, since M82-type sources will push into the sample 
and will not show up as luminous X-ray sources. In Fig- 
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urc 20, we again plot redshift versus total FIR luminos- 
ity, this time for the radio sources with LIRG luminosities 
(4 X 10'*'' < Lfir < 4 X lO''^ ergs s~*; vertical lines). 
The dashed horizontal line shows an adopted upper red- 
shift cut-off of z = 0.7. The ratio of highly-obscured AGN 
candidates to high X-ray luminosity sources for the LIRG 
population with 0.2 < z < 0.7 comes out to be 5.4, which 
can be taken as an extreme upper limit. 
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Fig. 19. — Number densities in two redshift bins (0.5 < z < 1.25 
and 1.25 < z < 2) of radio sources with either spectroscopic or 
photometric redshifts and Lpji^ > 4 X 10^^ ergs s~^, divided ac- 
cording to X-ray luminosity: diamonds denote sources with either 
Lo.5-2 keV or L2-S keV > 10^^ ergs s~^, and squares denote sources 
with both soft and hard X-ray luminosities < 10*^ ergs s~^. Squares 
are slightly offset for clarity. Poissonian Icr uncertainties are based 
on the number of sources in each redshift interval. 
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Fig. 20. — Redshift vs. total FIR luminosity for the radio sample 
with LIRG luminosities. Solid squares and solid circles (the latter 
arc CS sources) denote spectroscopic redshifts, and open squares de- 
note photometric redshifts. Large, open circles denote sources that 
have either io.5-2 keV or L2-S keV ^ lO'*'^ ^rgs s~^. The dotted 
curve shows the 60 /iJy completeness limit of the radio sample. The 
two solid vertical lines show the minimum luminosities of a LIRG 
and a ULIRG. The dashed horizontal line shows an upper redshift 
cut-off of z = 0.7 to provide as uniform a selection of LIRGs as 
possible, while still retaining a statistically significant sample. 

Models of the XRB generally require a population of 
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near-Compton-thick sources {Nh 
ical estimates are, on average, for 2 — 3 times as many such 
sources as X-ray detected sources (e.g., Fabian & Wors- 
ley 2004; Treister & Urry 2005). However, such obscured 



sources may preferentially lie at luminosities lower than 
quasar luminosities (e.g., Ballantyne ct al. 200G). Thus, 
our upper limit of 1.9 for the ratio of quasar- luminosity 
but not X-ray-luminous sources to quasar-luminosity, X- 
ray-luminous sources may be consistent with these results. 

If we assume that a roughly factor of two ratio of 
highly-obscured AGN candidates to high X-ray luminos- 
ity sources applies to the quasar luminosity population 
and shows no redshift effect, then we estimate that the 
accreted supermassive black hole mass density determined 
by Barger et al. (2005) for broad-line AGNs (their Eq. 9) 
will increase to 3.6 x lO'' Mq Mpc~'^. Although this is 
still within 2a of the local supermassive black hole mass 
density found by Yu & Tremaine (2002) of (2.9 ± 0.5) x 
10^ Mq Mpc^'^ for h =- 0.7, it does not leave much 
room for the obscured accretion from the optically nar- 
row AGNs. Applying the same factor of two corrections to 
Equations 6 and 7 from Barger et al. (2005), which give the 
accreted supermassive black hole mass density for all spec- 
tral types assuming bolomctric corrections, respectively, of 
85 and 35 for the optically narrow AGNs, would increase 
these numbers to 12 x 10^ and 6.3 x 10^ Mq Mpc"^. This 
could be viewed as a consistency check that there are not 
that many luminous obscured AGNs. 

10. CONTRIBUTIONS TO THE X-RAY LIGHT 

We would now like to understand what fraction of the 
X-ray light is coming from the different X-ray and radio 
populations, especially from the radio-identified ULIRG 
population. To study this, we performed source-stacking 
analyses on the CDF-N 2 Ms exposure, restricting to a 
9.5' radius circle rather than the 10' radius circle used 
elsewhere in this paper, since a few of the radio sources in 
the 10' radius circle lie just outside the X-ray image. 

For each sample, we first measured the X-ray fluxes of 
all of the sources in each of four passbands, 0.5 — 1, 1 — 2, 
2 — 4, and 4 — 8 keV, using the images and exposure maps 
given in Alexander et al. (2003). For the X-ray sample, 
we measured the fluxes of all of the sources in the Alexan- 
der et al. (2003) catalog that lie in the 9.5' radius region. 
(Note that Worsley et al. 2005 found an increase in total 
resolved flux of about 2 — 5% in the > 2 keV bands when 
they applied the systematic flux correction determined by 
Bauer et al. (2004) to account for Eddington bias and some 
additional aperture/photometry effects that were not con- 
sidered in Alexander et al. 2003, but we have not applied 
that correction here.) 

To calculate the contribution to the extragalactic back- 
ground light (EBL) from a given sample, we summed the 
fluxes in the 9.5' radius circle and divided by the corre- 
sponding area. We did not exclude any sources from the 
summations, even if a source were not significantly de- 
tected in a particular band or had a measured negative 
flux in the band. This simple averaging is not optimal, 
since the errors in the X-ray fluxes are smaller at smaller 
off-axis angles, but it avoids weighting the smaller central 
region more highly. 

To estimate the uncertainties on our surface brightness 
measurements, we first generated a grid of 60,000 points in 
the 9.5' radius area that are separated from each other by 
more than 4" (the aperture diameter used to measure the 
counts in the X-ray sources) and are located more than 6" 
away from any known X-ray source. We next measured 
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the fluxes in each of our four X-ray passbands at these 
60,000 positions. This serves as our base file. Then, for 
each sample for which we want to determine the uncer- 
tainties, we randomly generate a number (which matches 
the number of sources in the sample) of integers between 
1 and 60,000 and use the fluxes corresponding to those 
positions in our base file. We do this 1000 times, which 
gives us 1000 fltix values for each of the four passbands. 
We then order those values from the most negative to the 
most positive and use the 68% confidence interval to de- 
termine our uncertainties. (Most of our uncertainties are 
much smaller than the data points.) 




0.1 1.0 10.0 

Energy (keV) 



Fig. 21. — X-ray surface brightness vs. energy measured within a 
9.5' radius circle for the Alexander et al. (2003) X-ray sample (solid 
diamonds), for the X-ray sample without radio counterparts (open 
diamonds), for the present radio sample (solid squares), for the 
radio sample without X-ray— luminous counterparts (these sources 
have Z/0.5— 2 kcV or -L2— 8 keV < lO'*'^ ergs s^^) (open triangles), 
and for the radio sample without X-ray counterparts at all (open 
squares). The uncertainties arc 68% confidence intervals, as de- 
scribed in the text. For comparison purposes, the diagonal solid 
lines show the 1.4 photon index of the XRB normalized to the high- 
est energy bins of the X-ray surface brightness measurements for 
the radio sample, for the non-X-ray-luminous radio sample, and 
for the X-ray-undetected radio sample. The horizontal solid lines 
show the 2 — 10 keV XRB measurement of Revnivtsev et al. (2005) 
converted into narrower energy bands by adopting a photon index 
of 1.4. The thickness of the lines represents the uncertainties on 
their measurement. 

In Figure 21, we show the EBL from the X-ray sources 
in the Alexander et al. (2003) catalog with solid diamonds. 
Coincidentally, they match very well the 2 — 10 keV XRB 
measurement made by Revnivtsev et al. (2005) from a re- 
analysis of the HEA01/A2 data (thick horizontal lines; 
the thickness denotes the uncertainty range on their mea- 
surement ). Note that we converted the Revnivtsev et al. 
(2005) 2- 10 keV measurement to 2 - 4 keV, 4-8 keV, 
and 8—10 keV by adopting a photon index of 1.4. We 
chose this determination to compare with since the A2 or 
Cosmic X-ray Experiment instrument was especially de- 
signed for accurate measurements of the XRB over a very 
wide sky solid angle (e.g., Marshall et al. 1980; Boldt 1987; 
Gruber et al. 1999), and the design of its detectors allows 
the internal instrumental background to be separated from 
the XRB with almost absolute accuracy (Rothschild et al. 
1979; Boldt 1987). 

However, for this work, we are not interested in ac- 
curately determining the resolved fraction of the XRB. 



Such an analysis would require integrating the known 
log N — log S distribution to take into ac;count the rare, 
bright sources that are not sampled in this deep pencil- 
beam survey, extrapolating the number counts to take into 
account the sources with fainter fluxes, and deciding which 
of the XRB measurements to compare with. Indeed, this 
was done by Moretti et al. (2003), and there is very good 
agreement between their result and the Revnivtsev et al. 
(2005) measurement. Rather, it is sufficient for us to mea- 
sure the X-ray stirface brightnesses from the known X- 
ray sources in the CDF-N and see that they are roughly 
consistent with the IIEA01/A2 XRB measurement. In 
other words, we are not concerned with uncertainties at 
the 10 — 20% level. Our aim is just to see roughly what 
fraction of the X-ray light comes from the various popula- 
tions. 

We begin by considering only the Alexander et al. (2003) 
X-ray sources that are not radio sources in the Richards 
(2000) catalog. In Figure 21, we show the computed X- 
ray surface brightnesses from these sources with open di- 
amonds. We find that they contribute 50 ± 0.7% of the 
X-ray light in the 4 — 8 keV band. 

Next, we examine how much the total radio sample con- 
tributes to the X-ray light. In Figure 21, we show our com- 
puted X-ray surface brightnesses for this sample with solid 
squares. We see that the radio sample also contributes 
50 ± 0.5% of the X-ray light in the 4 - 8 keV band. In 
other words, radio sources to the Richards (2000) 1.4 GHz 
flux limit of 40 fiJy are markers for half of the X-ray light 
at these high energies. This contribution is the same as the 
contribution from the X-ray sources that are not detected 
at these radio fluxes. 

However, most of the light contributed by the radio sam- 
ple comes from the X-ray-luminous radio sources. To 
show this, we next measure the X-ray surface bright- 
nesses of the sample of radio sources without X-ray- 
luminous counterparts (these sources have io.5-2 kcV or 
L2-S kcV < 10^^ ergs s^"'^; hereafter, non X-ray luminous) 
in the Alexander et al. (2003) catalog (open triangles), 
as well as of the sample of radio sources without X-ray 
counterparts at all (hereafter, X-ray-undetected) (open 
squares). 

We find the contributions to the 4 — 8 keV light from 
the non-X-ray-luminous population to be 2.3 ± 0.4% and 
from the X-ray-undetected population to be 1.2 ± 0.3%. 
Thus, although some of the radio sources are contributing 
additional X-ray light that is not already in the X-ray sam- 
ple, such contributions are very small. This immediately 
tells us that the current radio source population cannot 
account for the background light that has been suggested 
may be missing at these energies. Indeed, the percent- 
ages are about a factor of 5 lower than those predicted for 
Nff ^ 10^^ cm~^ sources at these energies in typical XRB 
synthesis models (e.g., Comastri, Gilli, & Hasinger 2006). 

However, the inverse statement that some of these non- 
X-ray luminous or X-ray undetected radio sources con- 
tain highly-obscured AGNs is true. We can see this by 
comparing the shapes of the X-ray surface brightness mea- 
surements for the various samples with the shape of the 
XRB. To make it easier to do the comparisons, in Fig- 
ure 21, we normalize the 1.4 photon index of the XRB 
to the surface brightnesses in the highest energy bin for 
all but the X-ray sample and the X-ray sample without 
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radio counterparts (to avoid cluttering up the figure too 
imich) . We can sec that the shape of the surface brightness 
measurements for the radio sample is slightly softer than 
that of the XRB, while the shape of the surface brightness 
measurements for the X-ray sample without radio coun- 
terparts is slightly harder than the XRB. However, the 
shapes of the surface brightness measurements for the non- 
X-ray-luminous radio sample and the X-ray-undetected 
radio sample both appear much harder than that of the 
XRB. This result suggests that these two populations are 
consistent with containing highly-obscured AGNs. 




1.0 

Energy (keV) 



10.0 




1.0 

Energy (keV) 



10.0 



Fig. 22. — X-ray surface brightness vs. energy measured within a 
9.5' radius circle for three redshift intervals: < z < 0.5 (squares), 
0.5 < z < 1.5 (diamonds), and 1.5 < z < 3 (triangles) (a) for the 
radio sample without X-ray-luminous counterparts (these sources 
have io.6-2 keV or ^-2-8 keV < 10*^ srgs s"'^) in the Alexander et 
al. (2003) catalog (solid symbols) and (b) for the radio sample with- 
out X-ray counterparts at all (open symbols). Small offsets have 
been applied in the a;-direction to the diamonds and to the trian- 
gles for clarity. The uncertainties are 68% confidence intervals, as 
described in the text. For comparison purposes, the diagonal solid 
lines show the 1.4 photon index of the XRB normalized to the high- 
est energy bins of the 0.5 < z < 1.5 and 1.5 < z < 3 samples. The 
horizontal solid lines show the 2 — 10 keV XRB measurement of 
Revnivtsev et al. (2005) converted into narrower energy bands by 
adopting a photon index of 1.4. The thickness of the lines represents 
the uncertainties on their measurement. 

One might wonder how the shapes of the X-ray sur- 
face brightness measurements for the non X-ray luminous 
radio sample and for the X-ray-undetected radio sample 
vary with redshift. We show this in Figures 22a and b 
for the redshift intervals ^; = — 0.5 (squares), 0.5 — 1.5 



(diamonds), and 1.5 — 3 (triangles). The solid symbols in 
Figure 22a denote the non X-ray luminous radio sample, 
and the open symbols in Figure 22b denote the X-ray- 
undetected radio sample. Small offsets have been applied 
in the .x-direction to the diamonds and to the triangles for 
clarity. The uncertainties are again the 68% confidence 
intervals. We normalize the highest energy bins of the 
z = 0.5 — 1.5 and ^ = 1.5 — 3 samples to the 1.4 photon in- 
dex of the XRB for comparison. In both figures, we can see 
that the lowest redshift interval sample (squares) is quite 
soft due to the contributions from star-forming galaxies, 
while the higher redshift interval samples are hard due to 
the obscured nature of the sources. 
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Fig. 23. — X-ray surface brightness vs. energy measured within 
a 9.5' radius circle for the z < 2 radio-identified ULIRG sample 
with X-ray-luminous (1/0.5-2 keV or L2-8 keV > 10*'^ srgs s~^) 
counterparts in the Alexander et al. (2003) catalog (solid squares), 
without X-ray luminous counterparts (open triangles), and with- 
out X-ray counterparts at all (open squares). For comparison pur- 
poses, the diagonal solid lines show the 1.4 photon index of the 
XRB normalized to the highest energy bins of the X-ray-luminous 
and non-X-ray-luminous samples. The horizontal solid lines show 
the 2 — 10 keV XRB measurement of Revnivtsev et al. (2005) con- 
verted into narrower energy bands by adopting a photon index of 
1.4. The thickness of the lines represents the uncertainties on their 
measurement. 

Since in §9 we focused on the z < 2 radio-identified 
ULIRG sample, we now repeat our source-stacking analy- 
ses using only these data. Our results are shown in Fig- 
ure 23. We use solid squares to denote the X-ray surface 
brightnesses measured for the z < 2 ULIRG sample with 
X-ray-luminous (Lo.5-2 keV or L2-8 keV > 10'*^ ergs s~^) 
counterparts in the Alexander et al. (2003) catalog. These 
sources contribute 25 ± 0.1% of the X-ray light in the 
4 — 8 keV band, which shows that sources with ULIRG 
radio powers are responsible for contributing a very sub- 
stantial fraction of the XRB. We use open triangles to de- 
note our measurements for the z <2 non-X-ray-luminous 
ULIRGs. These contribute only 0.5 ± 0.2% of the X-ray 
light in the 4 — 8 keV band. Finally, we use open squares to 
denote our measurements for the z < 2 ULIRGs with no X- 
ray detections at all in the Alexander et al. (2003) catalog. 
These contribute only 0.4±0.2% to the 4-8 keV light. The 
uncertainties are again the 68% confidence intervals. We 
normalize the highest energy bins of the X-ray-luminous 
ULIRGs and the non-X-ray-luminous ULIRGs to the 1.4 
photon index of the XRB for comparison. We see that 
the shape of the X-ray-luminous ULIRGs is quite well 
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matched to the XRB, while the shapes of the non- X-ray 
luminous ULIRGs and the X-ray undetected ULIRGs are 
considerably harder, rising steeply towards the higher en- 
ergies. Again, this is consistent with these samples being 
highly obscured. However, even if we interpret them as 
such, their contributions to the 4 — 8 keV light are very 
small, and they are unlikely to contribute substantially to 
the XRB at even higher energies. 

11. SUMMARY 

In this paper, we analyzed the nature and evolution of 
micro Jansky radio sources using a highly spectroscopically 

complete, deep VLA survey of the HDF-N region. We 
supplemented the spectroscopic identifications with pho- 
tometric redshifts measured from the rest-frame ultraviolet 
to MIR spectral energy distributions. Our results are as 
follows. 

• We found that the fraction of radio sources that can 
be optically spectroscopically identified is fairly indepen- 
dent of radio flux, with about 60 — 80% identified at all 
fluxes. We spectrally classified the galaxies into four spec- 
tral types: absorbers, star formers, Seyfert galaxies, and 
broad- line AGNs. 

• We did not confirm the existence of an X-ray-radio 
correlation for star-forming galaxies in either the hard or 
soft X-ray bands. Previous claims of such a correlation 
appear to be the result of selection effects. 

• We did not observe any correlation between 1.4 GHz 
flux and R magnitude, which helps to explain the unifor- 
mity of the spectroscopically identified fraction with radio 
flux. Previous claims of the existence of such a correlation 
were likely due to the much shallower optical observations 
that were used in those studies. 

• We also did not observe any correlation between 
1.4 GHz flux and redshift for the spectroscopically and 
photometrically identified sources. The redshift distribu- 
tion is consistent with being invariant with radio flux, 
though we may be missing the higher redshift tail of the 
redshift distribution function due to optical spectroscopic 
bias. 

• The observed invariance in the rest-frame colors and 
only small amount of change in the absolute magnitudes of 
the host galaxies suggests that radio observations sample 
the same type of host galaxies at all redshifts and that the 
host galaxy properties are not evolving much with redshift. 

• There is dramatic evolution in the radio powers of the 
host galaxies with increasing redshift, and similar optical 
galaxies are hosting more powerful radio galaxies at higher 
redshifts. A quantitative description the evolution of the 
1.4 GHz luminosity function can be found in Cowie et al. 
(2004a). 



• Assuming that the locally determined FIR-radio cor- 
relation holds at high redshifts, we estimated total FIR 
luminosities for the radio sources. We note that these 
will be overestimates in the case of radio-loud AGNs, but 
we conservatively did not try to remove the radio-loud 
AGNs from the sample, since we are only determining up- 
per limits. For the radio sources with total FIR luminosi- 
ties comparable to quasar-like bolometric luminosities, we 
obtained an upper limit of 1.9 for the ratio of candidate 
highly-obscured AGNs to X-ray-luminous sources. 

• We measured the X-ray surface brightnesses for vari- 
ous X-ray and radio populations within a 9.5' radius circle 
to determine their contributions to the X-ray light. For 
the known X-ray sources in the CDF-N, we found our 
measurements to be (coincidcntally) consistent with the 
HEA01/A2 XRB measurement made by Revnivtsev et al. 
(2005). We found that the Richards (2000) radio sample 
contributes half of the 4 — 8 keV light, which is also how 
much the X-ray sources without radio counterparts con- 
tribute. However, most of the light from the radio sample 
comes from the X-ray-luminous radio sources. The non- 
X-ray-luminous and X-ray-undetected subsamples con- 
tribute only 2.3 ± 0.4% and 1.2 ± 0.3%, respectively. This 
tells us that the current radio source population cannot 
account for the background light that one might think is 
missing at these energies. However, the shapes of the sur- 
face brightness measurements for these two samples both 
appear much harder than that of the XRB, which suggests 
that these two populations are consistent with contain- 
ing highly-obscured AGNs. For the z < 2 radio-identified 
ULIRG sample, we found that the X-ray-luminous sub- 
sample contributes a very substantial one-quarter of the 
4 — 8 keV light, while the non- X-ray- luminous and X-ray- 
undetected subsamples each contribute less than a percent. 
Again, the shapes of the latter two subsamples are consid- 
erably harder than the XRB and hence are consistent with 
containing highly-obscured AGNs. 
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Table 1 

Bi-LOBAL Source and 11 Sources Without Obvious Counterparts 



Radio (J2000.0) Radio Flux Redshift 

R.A. Decl. (/iJy) 



12 


37 


25.73 


62 


11 


28.50 


5960"* 


1.641^ 


12 


37 


02.29 


62 


23 


31.50 


78.7 


0.642^^ 


12 


35 


40.39 


62 


16 


23.80 


60.1 




12 


35 


41.00 


62 


18 


28.20 


89.0 




12 


35 


57.71 


62 


08 


08.41 


68.6 




12 


36 


08.24 


62 


15 


53.00 


59.3'' 




12 


36 


24.28 


62 


10 


17.00 


54.2'* 




12 


36 


43.88 


62 


05 


59.20 


62.5 




12 


36 


46.70 


62 


12 


26.49 


72.0 




12 


36 


51.72 


62 


05 


02.49 


57.6'* 




12 


36 


54.70 


62 


10 


39.60 


48.2 




12 


37 


23.05 


62 


05 


39.59 


78.5* 





Note. — Units of right ascension are hours, minutes, and sec- 
onds, and units of dechnation are degrees, arcminutes, and arc- 
seconds. 

''Also present in the Biggs & Ivison (2006) catalog. 

'^Bi-lobal source centered on a small red galaxy with optical 
coordinates (12:37:25.9, 62:11:29.0). 

■^Source appears to be associated with a pair of interacting 
galaxies a,t z = 0.642. However, because the association is 
uncertain, we assume this source is unidentified. 
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Table 2 

(5^^'i']^52§iMb)8kiFTS Due to Uncertain iDENTiFicAft^g" Flux 


Literature 


R.A. 


Decl. (^Jy) 


Redshift 



12 


35 


53.26 


62 


13 


37.70 


58.4 


2.098'^ 


12 


36 


08.24 


62 


15 


53.00 


59.3 


0.459'^ 


12 


36 


11.41 


62 


21 


49.79 


111.0 


0.294^^ 


12 


36 


21.27 


62 


17 


08.40 


148.0 


1.992'^ 


12 


36 


36.91 


62 


13 


20.41 


50.0 


0.680<^ 


12 


36 


42.10 


62 


13 


31.41 


467.0 


4.420'' 


12 


36 


46.70 


62 


12 


26.49 


72.0 


2.970^ 


12 


36 


51.76 


62 


12 


21.30 


49.3 


0.401S 


12 


36 


56.60 


62 


12 


07.60 


46.2 


0.321=^ 


12 


37 


01.57 


62 


11 


46.60 


128.0 


0.884'^ 


12 


37 


06.77 


62 


07 


22.50 


72.8 


0.518<i 


12 


37 


50.27 


62 


13 


59.00 


90.5 


0.231<i 



Note. — Units of R.A. are hours, minutes, and seconds, and units of Decl. are degrees, arcminutes, and arcseconds. 

'^Chapman et al. (2005) redshift for SMM J123553.26 + 621337.7 and Swinbank et al. (2004) redshift for SMM 
J123621.27 + 621708.4. Each radio source has an 8 /ini counterpart, but there is no blue hght at those positions. The 
redshifts appear to be of neighbor objects, which may or may not be associated with the radio sources. 

''Wirth et al. (2004) redshifts for GOODS J123608.02 + 621554.0 and GOODS J123701.80 + 621144.2. The first radio 
source is far off-axis. The second radio source has an 8 /xm counterpart, but the redshift is of a neighbor galaxy. 

^Barger et al. (2002) redshift for a neighbor galaxy, CXOHDFN J123611.40 + 622149.9. 

^Barger et al. (2000) redshifts for neighbor galaxies. 

**Waddington et al. (1999) redshift. The radio source has an 8 fim counterpart, but the identification is offset. 
^Steidel et al. (2003) redshift. The radio source has a weak 8 fim counterpart, but the redshift is of a neighbor galaxy. 
^Lanzetta et al. (1996) redshift. The radio source has an 8 /^m counterpart, but the redshift is of a bright neighbor. 



